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Abstract 

We consider the totally asymmetric exclusion process in discrete time with generalized updating 

rules. We introduce a control parameter into the interaction between particles. Two particular 

values of the parameter correspond to known parallel and sequential updates. In the whole range 

, . of its values the interaction varies from repulsive to attractive. In the latter case the particle flow 

■^ ' demonstrates an apparent jamming tendency not typical for the known updates. We solve the 

^' 

' master equation for N particles on the infinite lattice by the Bethe ansatz. The non-stationary 

_L ' solution for arbitrary initial conditions is obtained in a closed determinant form. 
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I. INTRODUCTION 

The totally asymmetric simple exclusion process (TASEP) is a paradigmatic model of 
stochastic systems of interacting particles demonstrating non-equilibrium behavior l|-|3|. 
The discrete-time dynamics of the model on the integer one-dimensional lattice is charac- 
terized by one of possible updating rules. The most important cases are the backward- 
sequential, parallel, sublattice-parallel and random sequential updates [4]. The first two 
of them are investigated in more detail than others. For a finite number of particles, the 
dynamics can be defined through a master equation 

P(X,t + l) = ^T(X,X')P(X',t) (1) 

X' 

where components of A^- dimensional vector X = {xj}, (xi < X2 < ■ ■ ■ < xn) are positions of 
A^ particles and T(X, X ) is the transition probability to go in one time step from configura- 
tion X to configuration X. For the backward-sequential update, each particle may take one 
step to the right with probability p if the target site is vacant at the beginning of the time 
step or becomes vacant at the end of the time step due to motion of the particle in front. 
For the parallel update, the motion to the right is allowed only if the target site is vacant 
at the beginning of the time step. By iterating ([T]) one obtains the solution of the master 
equation for any given initial configuration X", that is the conditional probability to find 
a particle configuration X at time step t, given that the process started from configuration 

A common property of the above mentioned updating rules is the short range repulsion 
between particles due to the exclusion condition. On large scales, the particle density evolves 
according to the Burgers equation [5I|. Fluctuations of the density and current on an appro- 
priate scale are the subject of intensive studies [6|-[8| over the past two decades. The second 
important property of the mentioned updates is their "solvability" which means that the 
master equation ([1]) can be solved by the Bethe ansatz method and the Green functions can 
be calculated explicitly [9|, [lO|. In this paper, we propose a generalization of updating rules 
in order to modify the first property retaining the second one. Namely, besides the exclusion 
condition, we allow an attractive interaction between neighboring particles which changes 
evolution of the system drastically. 

Specifically, we introduce a control parameter into the interaction between particles. 
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FIG. 1: The probabilities of different jumps of two neighboring particles. 

Depending on its value, we obtain an attraction or an additional repulsion between particles. 
We consider the TASEP on the infinite one- dimensional lattice and solve the master equation 
for arbitrary initial conditions to find the Green function in the determinant form as it was 
done in Ill4l3l|. The solution generalizes the known results for the backward-sequential and 
parallel updates and coincides with them for particular values of the parameter of interaction. 
The article is organized as follows. In the Section II, we formulate the model and define 
basic notions needed for the solution. In the Section III, we explain a peculiarity of the 
master equation in our case and transform the evolution operator to a form appropriate for 
the Bethe ansatz. Solving the master equation for arbitrary initial conditions, we obtain 
the closed expression for the Green function. The Appendix contains the proof of the Bethe 
ansatz in the general A^-particle case. 

II. DEFINITIONS 

Consider A^ particles moving to the right on the one- dimensional lattice. At each discrete 
moment of time, each lattice site can be occupied by at most one particle. To describe 
the dynamics of the model, we follow the backward-sequential rule: at given time step, 
we scan the particle configuration from the right to the left and consider situations when 
two particles meet in neighboring sites just before the time step (see Fig{T]). If the right 
particle of the pair does not hop at given step, the left particle stays at its position with 
probability 1. If the right particle hops, the next particle jumps to the right with probability 
(1 + z/)p, or stands with probability 1 — (1 + v)p. A natural restriction on the parameter v is 
< (1 + z/) < 1/p. For z/ = 0, we have the usual TASEP with backward-sequential update. 
For V = —1, we obtain the TASEP with parallel update. 

If 1/ > 0, the system demonstrates a collective behavior which is quite typical for the car 
traffic. To see this, consider a limiting illustrative example z/ = l,p = 1/2. Isolated particles 
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FIG. 2: Typical behaviour of particle space-time trajectories at (a) rarefaction regime, withp = 0.1, 
ly = —1, and (b) jamming regime, p = 0.1 and ly = 8.95. The dynamics of A^ = 100 particles 
is simulated for T = 2000 time steps starting with the initial configuration where particles are 
alternating with empty sites. 

move forward as usual with the probabihty p = 1/2. If two particles turn out nearest 
neighbors, they become "stuck together" and move synchronously because the situation 
when the right particle moves and the left one stands is forbidden. This micro-jam grows 
involving into oneself new and new particles until the moment when all particles move 
synchronously in a huge A^-particle jam. Thus, when varying the parameter u from —1 to 
{1/p— 1) the system passes through a qualitative change of particle density behaviour: from 
rarefaction to clumping. Two typical examples of such behaviour are shown in Fig. ([2]) 

We start solving the problem with several definitions. We introduce the Hilbert space 
supplied with the complete left and right bases consisting of vectors (X| and |X) respectively, 
where X runs over all particle configurations, with inner product 



{X\X') = S{X,X'). 



(2) 



The state of the system at any time step can be associated with the state vector 



\P,) = J2PtiX)\X). 

{X} 

In terms of the state vectors, the master equation ([1]) takes a simple operator form 



(3) 



l^m) = T|P,) 



(4) 



where the evolution operator T is defined by the transition probabihties 

T= Y. \X)T{X,X'){X'\. (5) 

The conditional probability P{X,t\X^,0) can be represented as the matrix element 

P(X,t|X°,0) = (X|T*|X°). (6) 

To evaluate matrix elements, we construct the set of left eigenvectors \Bz) of the operator 
T 

T\Bz) = A{Z)\Bz) (7) 

and the adjoint set of right eigenvectors {Bz\ 

{Bz\T = A{Z){Bz\, (8) 

where Z is a A^- dimensional parameter. Now, the transition probability ([6]) can be reduced 
to the evaluation of the sum 

{X\ T* |X°> = Y, {X\T'\Bz) {Bz\X') = J] A* (Z) {X\Bz) {Bz\X') , (9) 

z z 

where the range of summation is to be defined from the boundary conditions. 

III. SOLUTION 

To understand a peculiarity of the eigenproblem ([7]) in our case, we consider first the 
two-particle problem. If two particles are at least two empty sites apart from each other, 
they move independently and Eq.Q is equivalent to the equation 

2 

A{Z){x,,X2\Bz) = Y. Up'' ^^ - P^'~'' 

{kuk2}i=l 

X (xi - ki,X2- k2\Bz) . (10) 

where the summation is over A;i, k2 which take values 1 or depending on whether the particle 
decided to jump or not. If particles are neighbors, xi + 1 = X2, the terms in the RHS of the 
actual (interacting) equation ([71) are p^(l + v) {xi — 1, 0:2 — l\Bz), p{l — p) {xi — l,X2\Bz), 
(l—p) {xi,X2\Bz), whereas the termp(l— p) (a;i,X2 — l|-Bz) = p(l— p) {xi,xi\Bz) in the free 



equation (fTOj) is beyond the set of allowed configurations. The strategy of the Bethe ansatz 
is to define formally the forbidden term via allowed terms in such a way that Eq.Q retains 
its "free" form f IlOp in the case Xi + 1 = X2. Then the formal relation between forbidden 
and allowed terms can be used for determination of coefficients of the Bethe substitution 

(a:i,X2|5z) = ^12^1""^ ^2""' + A^.zl^^^z^^^ (11) 

This strategy is effective for the TASEP with the backward-sequential update. In our 
case (as well as in the case of parallel update [l2|) an additional problem arises. To see 
it, consider configuration {xi,X2} where xi + 2 = X2- The terms in the RHS of Eq.([7]) 
are p^ {xi - 1,X2 - l|i?z), p(l - p) {xi - l,X2\Bz), (1 - p^ {xi,X2\Bz) and p{l - (1 + 
i')p) (xi,X2 — l\Bz). The last term differs from its free version p{l — p) (xi,X2 — ^\Bz) 
and we have no forbidden terms like p(l — p) (xi, xi\Bz) in the previous paragraph, to com- 
pensate this difference. Instead, we can notice that the configuration {xi,a;2 — 1} contains 
the pair of neighboring particles, which disappears in the configuration {xi,a;2}. Thus, if 
we multiply the transition probability T(X, X ) by the factor A = (1 — p)/(l — (1 + v)p) 
corresponding to dissociation of the pair, we restore the free form of the last term in the 
RHS of fITOl) . In the A^-particle case, the factor A should be ascribed to each dissociated 
pair. 

More formally, we define an auxiliary matrix 

Tn = VTV-^ 



where D is a diagonal operator 



^_.^ |X)(X| 
{X} ^ ' 



and iy(X) depends on the number of pairs Np{X) in configuration X: W{X) = \'^p^^>. 
Operators T and To have the same eigenvalues A^ and the right eigenvectors are related 

by 

\Bz)=V-'\Bl). (12) 

The matrix elements of Tq are 

To(X,X') = ^^T(X,X'). (13) 



The two-particle eigenproblem for operator To has the free form 

K{z){x,MB'z) = Yl Up'' (^ - P)"'"' 

X (xi - ki, X2 - k2\B'^z) ■ (14) 

if the distance between particles exceeds 1. If particles are neighbors, a;i + 1 = X2, the 
terms in the RHS of ([H]) are p'^{l + v) (xi - 1,X2 - l\B%), p{l - p)/X {xi - 1,X2|5|), 
(1 — p) {xi,X2\B^) and the term p{l — p) (a;i,a;2 — 1|-B^) in the free equation is forbidden. 
Imposing the condition 

{x,x\B°z) = ^^i{x-l,x\B°z)- (x-l,x + l|fi|)) + (x,x + l|fi|). (15) 

1 — p 

on the forbidden term, we convert the remaining three terms into the free form, ensuring 

correctness of flT^ for arbitrary allowed configurations {xi,X2}- 

In Appendix, we prove that (ITSl) written for the general case as 

(...,a:,x,...|50) = ^((...,x-l,x,...|fii)-(...,a;-l,x + l,...|i?|)) 

+ (..., x,x + l,... 15°) 

is sufficient to convert the general A^-particle problem into the free form 

N 

A{Z){x„...,x^\B'z) = EnP'"(l-P)'^'' 

{k,} i=i 

[xi- ki,...,XN - kN\Bz) . (17) 



X 



After that we can look for B^ (xi, . . . , xn) in the form of the Bethe ansatz, 

(xi, . . . , xm\B'z) = Y, A^.-^^.^T ■ ■ ■ C"' (18) 

where parameter Z is a set of A^ complex numbers Z = {zi, . . . , zn}- The summation is 
over all permutations a of numbers 1, . . . , A^. The substitution of ( TT8|) to the equation ( ITTll 
leads to the eigenvalues 

N 

A{Z) = l[{l-p + pz,), (19) 

and the substitution of (fT8l) into the constraint flT6l) gives the relation between the amplitudes 
^cri...o-jv5 which differ from each other only by permuted indices 

A...ij... _ 1 - l/zi 1 - 72;^- 



A...ji... 1 - l/zj 1 - -fZi 
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(20) 
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FIG. 3: The transition probability for the processes defined by operators T and Tq. 



where 7 



1/A. Separately, amplitudes Ao-i...o-jv ^-^e 

AT 



A. 



^^{W) 



o-i...crjv 



n 



4 = 1 



7^<7, 



1/^., 



(21) 



where V ({o"}) is the parity of a permutation {a}. Then, we are in position to write the 

element {X\Bz) in ^. 

To find ^i?^|X°), we examine the right eigenproblem ([8]). As above, we start with the 

two-particle case. If X2 — a;i > 2, Eq.(l8]) is equivalent to the equation 

2 



X {Bz\xi + ki,X2 + k2) . 



(22) 



where the summation is over fci, ^2 taking values 1 or 0. If X2 — a^i = 1, the RHS of ([8]) 
contains three terms: p^(l + z/) (^Bz\xi + 1,X2 + l), p(l — (1 + z/)p) (_B^|xi,X2 + l) and 
(1 — p) (^Bz\xi,X2). At the same time, the RHS of Eq. (l22i) contains four terms in this case, 
and one of them, p{l — p) (^Bz\xi + 1^X2), is forbidden. We can see from Figj3] 

that the right eigenproblem for the operator T in the two-particle case is similar, up to 
the inversion of the coordinates, to the left eigenproblem of the operator Tq. This is true for 
the general A^-particle case and, therefore, we can use for (^Bz\X^ the similar Bethe ansatz: 



XI XN 

^o-i...crjv^CTi • • • ^o-jv 



(23) 



W} 



with the amphtudes A„ given by 
Collecting obtained expressions, we have 



{X\Bz) = W^(X)detB, (24) 

{Bz\X) = ^detB (25) 

where the matrix elements Bij and Bij of the matrices B and B are 

The factor 1/A^! is a normalization constant which follows from the condition 

J2{X\Bz){Bz\X) = l. (27) 

z 

The range of summation over Z is defined from the boundary conditions. For the case of 
the infinite lattice treated here, the spectrum is continuous, such that the sum in (Q can be 
written as integral 

where each Zi runs along the closed contour Tqi in complex plane, which encircles the points 
2; = 0, 1 and leaves all other singularities outside. The choice of the contour is justified by 
the completeness requirement in a sense that the set of operators {\Bz) {Bz\ '■ Z G r^;^} 
provides the resolution of the identity operator 

Y,{X\Bz){Bz\X') = {X\X'), (29) 

z 

which has been proved in [12] for the case of parallel update. The proof in the general case 
differs by technical details only. 

An explicit expression for the conditional probability P(X, t\X^, 0) follows from sequential 
substitutions ([26]) into ([25]), then ([25]) and ([19]) into ([9]) and ([6]). 

N 

(X|T*|X°> = iy(X)5^(-l)^(M)-Q^_^(^^_^o^^^)^ (30) 

{a} i=l 



where the function / is 

Finally, we obtain the Green function as the determinant of matrix M: 

P(X, t|X°, 0) = A^''^^) det M, (32) 

with matrix elements 

Mij = Fi_j{xi-x],t). (33) 

In the following part of this work, we will use the obtained Green function for calculation 
of correlation functions and probability distribution of the position of the last particle which 



is the Tracy- Widom distribution (ij, |l5| in the case of standard updates 
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Appendix A: induction 

In this section, using an induction we prove that using f lT6|) one reduces the problem with 
interaction to the free equation ( !T7|) . In fact below we give the proof only for the particular 
case of a single isolated cluster of particles. When there are several clusters which are more 
than one site apart from each other, they are independent, and this proof is applied to each 
clusters separately. In principle one has to consider also the case of clusters separated by 
single empty site. However technically this proof is similar to the one below and therefore 
we omit it. 

We start with the eigenproblem flT7|) for three particles situated in the cluster xi = x,X2 = 
X + 1,X3 = X + 2: 

3 

{k,} i=i 

x{xi-ki,X2-k2,X3-h\B'^z) ■ (Al) 
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We have to prove that condition (fT6!) converts (lAip into the eigenproblem of operator Tq. 
The RHS of flAl|) contains four forbidden terms {x,x + l,x + 1\B^), {x — l,x + l,x + 1|-B^), 
{x,x,x + 2\B^) and {x,x,x + 1|-B|). Using a straightforward extension of (IT^ . we define 
two of them as 

{x,x,x+2\B^z) = ■T^^{{x-l,x,x+2\Bl)-{x-l,x+l,x+2\Bl)) + {x,x+l,x+2\Bl) (A2) 



(x-l,x+l,x+l|S°) = — ^((x-l,x,x+l|S^)-(x-l,x,x+2|S^)) + (x-l,x+l,x+2|5|) 



For (x, X + 1, x + 1|-B|) and (x, x, x + Ij-B^-) we have the hnear system 

(x,x + l,x + l|E|) = j^((x,x,x+l|5^) -(x,x,x + 2|5|)) + (x,x + l,x + 2|S|) 



X, 



x,x + l|5|) = Y^((x-l,x,x + l|fi|) - (x-l,x + l,x + l|50)) + (x,x + l,x+l|fi|) 



which gives 

(x,x + l,x + l|5°) = , ^ ^ , r(x-l,x,x + l|E|)-(x-l,x,x + 2|M)) + (x,x + l,x + 2|E|) 

(A4) 
and 

(x,x,x+l|5|) = -^((x-l,x,x+l|5|)-(x-l,x+l,x+2|5|)) + (x,x+l,x+2|5|) (A5) 
1 — p 

Using definitions (jM]) , dMl) , dAll) and (jMl) we rewrite ([SI]) as 

Az{X\B°z) = (1 + «^) V(a; - 1, X, X + l\B°z) + Qi^^ x + 1, x + 2|50 ) + 

(1 + y)p^{q — yp){x — l,x, X + 2\B^z) + p{q — '^p){x — l,x + 1,x + 2\B%) 

(A6) 

what corresponds exactly to the eigenproblem of operator Tq for the cluster X = {x, x + 

l,x + 2}. 

Using an induction, we can generalize flA4p and (jASj: 



x,x + l,x + 2, ...,x + A^-2,x + iV-2|5|) = (x,x + 1, x + 2, ..., x + A^ - 1|5|)- 

(A7) 



^N-lpN-l 



(i-p)^-i ((a: - 1, X, X + 1, ■■■, X + A^ - 2|g|) - (x - 1, x,x + 1, ..., x + A^ - 3,x + A^ - 1|5|)), 



(x,x,x + l,...,x + A^-2|5|) = (x,x + l,x + 2,...,x + A^- 1|S|) + 

/ \ (A8) 

Y^f(x-l,x,x + l,...,x + A^-2|B|) - (x-l,x + l,x + 2,...,x + A^-l|5|)j. 
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Suppose now that for A^ particles we have 

11 1 

fci=0fe2=0 fcjv=0 

{x - ki,x + I - k2, ...,x + N - I - kN\Bl) = {I + u^-^p^^x -l,x,...,x + N - 2\B%) + 

N-l 

{I - p){x,x + I, ...,x + N - l\B%) + E {l + yf~^p^{q-yp)^ 

k=l 

(x — l,x,x + l,...,x + A; — 2,x + A;,x + fc + l,...,x + A^ — 1|-B|) 

(A9) 

Then, for A^ + 1 particles we obtain 

11 1 

A:i=0/C2=0 fcjv+i=0 

(x - fci,x + 1 - /c2, ■■■■,x + N - kN+i\B%) = 
1 
J2 pkN+i(^i _ py-kM+1 (^(^i ^ ,j)N-ipN X (^^_i^^^___^^^j^_2^x + N-kN+i\Bl) + 

A:jv+i=0 

N 

{l-p){x,x + l,...,x + N - l,x + N -kN+i\B9.) + Y.{1 + u)''-^p''{q - up)x 

fc=i 

(x- l,x,x + l,...,x + k,x + k + 2,x + k + 3,...,x + N - l,x + N - kN+i\B^z)). 

(AlO) 

Applying (lA7p .( lA8|l we derive 

(1 + u)^-^p^ [p{x -l,x,...,x + N -2,x + N - 1|50) + 

{l-p){x-l,x,...,x + N -2,x + N\B'^z)] + 
(1 -p)[il -p){x,x + l,...,x + N - l,x + N\B^z) +p{x,x + l,x + 2,...,x + N\B'^z) + 

^J^^((x - 1, X, X + 1, ..., x + N- l|fiO) _(^x-l,x,x + l, ..., X + A^ - 2, X + Ar|5| ))] + 



N 

E (1 + yY^^p^{q - vp) [(1 - p) X 
fc=i 

(x- l,x,x + l,...,x + A;-2,x + A;,x + fc + l,...,x + A^- l,x + A^|E|) + 

(Iz^[^^"V^+'((^ - 1, x, X + 1, ..., X + AT - 115° ) + 

Y^((l - p)^+i-V(a; - 1, X, ..., X + A; - 2, X + fc, X + A; + 1..., X + Ar|E|))- 

z/^-^p^(l - p)(x - 1, X, X + 1, ..., X + A^ - 2, X + A^|5|))]] 

(All) 

The coefficient of the term (x — l,x, ...,x + A^ — 2,x + A^ — Ij-B^') i^ 

N 

^ '^^ fc=i ^ *^^ (A12) 
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The coefficient of the term {x,x + 1, ...,x + N — l,x + A^|-B|) is 

a-pr+a-p)j^^Y^i^-pf^' = ^-p, (A13) 

The coefficient of the term {x — l,x, ...,x + N — 2,x + N\B^) is 

(1 + ur-^p^{i-p) - ^^-p""^' E (1 + ^r-\<i - ^p)(Tz^^^-^ = 

and the coefficient of the term {x — l,x, ...,x + k — 2,x + k,x + k + 1, ...,x + iV|_B^) is 

^ ^ ^ -(g - -P)/(1 + T^) = n_tw+i-J ^ - -^)^' (A15) 



Collecting all terms, we obtain for A^ + 1 particles 



1 1 



fci=0fc2=0 A;jv+i=0 

{x — ki,x + 1 — k2, ...yX + N — kN\Bz) = 
(1 + z/)>^-i(x - 1, x, ..., a; + iV - 1|50) + (1 - p)(x, a; + 1, ..., a; + Ar|50) + 

N 

E(l + /^)^"V^(g - up){x -l,x,x + 1, ...,x + k -2,x + k,x + k + 1, ...,x + N\B^), 

k=l 

(A16) 



q.e.d. 
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